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Abstract 


An iron beam experiment recently conducted at the Lawrence 
Berkeley Laboratory Bevalac by Benton et al. provides an opportunity 
for verifying the new Green’s function computer code (GRNTRN) and 
assessing the related nuclear database . The iron beam with 600 MeV/A 
at extraction traversed a series of beam transport elements, lead foil , 
and several triggering devices before impacting the target. Of these, 
only the 2.24 g/cm 2 lead foil and target are considered in the trans- 
port analysis with an assumed (inferred from calibration) beam energy 
of 557 MeV/A. A thin layer of CR-39 plastic nuclear track detector 
(PNTD) was placed in front of the target to monitor the incident flux 
and a stack of four PNTD’s placed behind to measure the linear energy 
transfer (LET) distribution of the transported beam. Test data are ana- 
lyzed for three separate targets: 2 g/cm 2 aluminum , 5-crn polyethylene , 
and 8-crn polyethylene. The two-layer GRNTRN results were mapped 
into the detector response function for comparison with the measured 
LET spectra. Reasonable agreement is obtained. Future research and 
analysis can be improved by using a more accurate isotope set or in- 
cluding other important media which significantly alter the beam. The 
assumption that the fragmentation cross sections are too small for alu- 
minum is inconclusive because the production of fragments by 2 g/cm 2 
aluminum can be more significantly affected by the fragment contribu- 
tion originating from the media in front of the target than that currently 
considered in the analysis. 

Introduction 

In designing a spacecraft for piloted missions or 
a commercial high- altitude transport aircraft, con- 
sideration must be given to protecting the crews 
(refs. 13) and passengers (refs. 4 10) from expo- 
sure to harmful radiation originating in space. This 
safety consideration requires development of highly 
efficient shielding computer codes which are practi- 
cal for integrated system design; a set of shielding 
computer codes (refs. 11 13) is being developed at 
the Langley Research Center. Most recently, a new 
method of using Green’s function for the solution 
of the heavy ion transport equations has resulted in 
a code that is not only efficient for engineering de- 
sign application but also suitable for a monoenergetic 
beam source which could be validated in the labo- 
ratory. (See refs. 14-17.) Because of the dearth of 
experimental measurements of fragmentation cross 
sections for heavy ion collisions, much uncertainty 
is associated with the present nuclear cross section 
database in the computer codes. This uncertainty 
has had a significant impact on the study of shield- 
ing. (See ref. 18.) 

An iron beam experiment, recently conducted 
at the Lawrence Berkeley Laboratory Bevalac by 
E. V. Benton using CR-39 plastic nuclear track 


detectors (PNTD’s) offers an opportunity for verify- 
ing the new Green’s function computer code (refs. 16 
and 17) and assessing the existing nuclear database. 
The linear energy transfer (LET) spectra of the at- 
tenuated iron beam and ion fragments were obtained 
by the PNTD’s which were placed behind an alu- 
minum or polyethylene degrader. Theoretical re- 
sults obtained by using the Green’s function com- 
puter code were compared with the measured LET 
spectra. 

Experiment 

Iron beams were accelerated to a 600-MeV/A 
nominal extraction energy. After extraction, the 
beam traversed a series of beam transport elements, 
lead-scattering foils (i.e., beam spreader), and several 
triggering devices before impacting the target. A se- 
ries of tests were performed by Benton et al. (Univer- 
sity of San Francisco) using various thicknesses of the 
target materials shown in table I; only the data analy- 
sis for the 2 g/cm 2 aluminum, 5-cm polyethylene, and 
8-crn polyethylene targets were completed for this 
study. Figure 1 shows the experimental arrangement 
of the lead foil, target, and CR-39 PNTD’s consid- 
ered in the current transport analysis. The remaining 
media in front of the target and PNTD’s which will 



affect energy and composition of the beam are being 
neglected because the 0.198-cm-thick (2.24 g/cm 2 ) 
lead foil is the major factor in changing the beam. 
The 557-MeV/A beam energy used for current anal- 
ysis is inferred from calibration tests by assuming 
the simple layout in which only the lead foil effect 
has been considered. The energy spread of the beam 
has a standard deviation of 0.2 percent. 
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subject to the boundary condition 


A thin layer of CR-39 PNTD placed in front of 
the target monitors the incident beam intensity and 
a stack of four CR-39 PNTD’s behind the target mea- 
sure the transported LET distribution. Calibration 
of these detectors with various ion beams of known 
LET resulted in a detector response function which 
is approximately Gaussian with an LET-dependent 
r (full-width at half-maximum) as shown in figure 2. 
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The above equations can be simplified by trans- 
forming the energy into the residual range as 
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Green’s Function Methods 


Transport Equations and Conventional 
Approach 

The transport equation for high-energy heavy 
ions is usually simplified by assuming the straight- 
ahead approximation and neglecting the target sec- 
ondary fragments (ref. 19) and is written as 


dx 
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where <fij(x,E) is the ion flux at x with energy E in 
MeV/A } Sj(E) is the change in E per unit distance, 
oj is the total macroscopic absorption cross section, 
and Ojh is the macroscopic cross section for the 
i collision of an ion of type k to produce an ion of 
type j. The solution to equation (1) is subject to the 
boundary condition 


0j(O, E) — fj{E) (2) 

which for laboratory beams has only one value of j 
for which fj(E) / 0 and that fj(E) is described by 
a mean energy and an energy spread a such that 
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By defining new field functions as 
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and with the solution of the ion fields given by 
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Note that v 3 is the range scale factor as i sjrj — v m Tm 
and is defined as Vj = Zj/Aj. The solution to equa- 
tion (11) is written as a perturbation series 
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The solution to equation (1) is given by superposition 
of Green’s function Gj ^ as 
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where ^(*,^,0 = 0 unless 

— (rj +x)<r' m <-j-Tj+x (17) 

Vm ^ rn 

for Vm>Vy If Vj>Vm, which can happen in neutron 
removal, the negative of equation (16) is used and the 
upper and lower limits of equation (17) are switched. 
The higher order terms are approximated as 


Reference 17 extends the method to the heavy ion 
transport in multilayered materials. In this new ap- 
proach, Green's function is constructed from a convo- 
lution product of very thin shield solutions whereby 
the higher order terms can be neglected. 

First, recall that the g function of n arguments 
was generated by the perturbation solution of the 
transport equation neglecting ionization energy loss 
(ref. 20) given by 
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Note that the approximation of r' Tn ) by 

equation (IB) is purely dependent on x and for i > 0 
which is represented as Q^ n {x). From above and 
reference 3, the solution to equation (1) is 

tpj(x, rj) = exp {-crjx)}j{rj + x) 
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where r' mu and r\ n( are given by the upper and lower 
limits of the inequality (17) and Fm{ r m) re f crs f° If 10 
integral spectrum 
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New Approach With Nonperturbative 
Method 

The higher order terms (third and above) of the 
perturbative solution derived from equation (21) may 
be important in the practical applications which re- 
quire calculations for deep penetration. Because 
computation of the higher order terms is inefficient 
particularly for the fragmentation of heavier projec- 
tiles, the perturbative approach given above is not 
suitable for engineering design applications. The 
following describes a new nonperturbative approach 
developed recently by Wilson et al. (See ref. 16.) 
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for any positive values of x and y. Equation (26) 
may be used to propagate the function gjm(x) over 
the solution space from very thin shield solutions. 
Equation (14) is then rewritten as 

Sjmix, rj, r' m ) « cxp(-crjx)6j m 6(x + r, - r' m ) 
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and the approximate solution of equation (1) is given 

by 
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Green’s Function Methods in a Shielded 

Medium 

The major simplification of the Green's function 
method results from the fact that the scaled (in trans- 
formed variables) spectral distribution of secondary 
ions to a first approximation depends only on the 
depth of penetration predicted in equations (16), 
(18), and (27). The first approach to a multi- 
layered Green’s function will rely on this observation 
and assume its validity for multilayered shields. 
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Consider a domain labeled as 1 which is shielded 
by another labeled as 2; the number of ions of type j 
at depth x in domain 1 due to ions of type m incident 
on domain 2 of thickness y is 

9i2j,n (*> y) = £ 9 hk (X) g 2km (y) (29) 
k 


The leading term in equation (29) is the penetrat- 
ing primaries and all higher order terms are within 
the bracket of 

9\2jm( x ^y) ~ cxp(— a \jx - a 2j y)S jm 

+ [9i 2 jm( x ^y) - exp {-<Tijx - (T 2j y)6 jm ] 

(30) 

The first term of the scaled Green’s function is 
then 

g i2j7n ( x ' V’ r J ’ r m ) 

= exp(-a lj x - (T 2 jy)6j m 6[x + rj - (r' m - py)} 

, . (31) 

where p is the range scale factor for the two media 
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A single value is assumed for p corresponding to 
600 MeV/A. The secondary contribution is similarly 
found by noting that equation (17) becomes 


~( r j +x+py)< 4 < rj + x + py (33) 

rn 

from which the average spectrum is evaluated. The 
complete approximation of Green’s function is then 


Gl2jm( T -y> r j’ r ,n) ~ exp(-cq jX - V2jy)t>jmH x + py + Tj - r' m ) 

+ v ± [g|2j m V ) ~ e *P( -°\jx - <r 2 j 

(x + py)(v m - i/j) 

(34) 

Equation (34) is our first approximation of Green’s 
function in a shielded medium of two layers and is 
easily modified to multiple layers. 

For the first spectral modification, the first colli- 
sion term has the properties 
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These properties are used to correct the average 
spectrum to ^ 
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is the midpoint of r' m between its limits given by 
equation (33). The b jm term of equation (36) has 
the property that 



(38) 


which ensures that the first term of equation (36) 
is indeed the average spectrum as required. The 
spectral slope parameter is 


b jnA x >y) 
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A similarly simple spectral correction could be made 
to the higher order terms. The spectral correction 
given in equation (39) is included in the present 
Green’s function computer code. 

Solution for Laboratory Beams 

The boundary condition appropriate for labora- 
tory beams is given by equation (3). The cumulative 
spectrum is given by 


F,(E) = i 
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The cumulative energy moment needed to evaluate 
the spectral correction is 
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The average energy on any subinterval (£), £ 2 ) is 
then 
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The beam-generated flux is 

rj) = exp(-(7]jX - o<i jy)fj(' r j 4- x + py) 
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where E is evaluated using equation (42) with E\ 
and E 2 as the lower and upper limits, respectively, 
associated with r f mi and r f mu . 

Nuclear Data Base 

The nuclear absorption and fragmentation cross 
sections needed for the transport calculation are 
generated by the Langley nuclear fragmentation 
(NUCFRAG) computer code (ref. 21) by using a re- 
duced set of 80 isotopes as listed in table II. In the 
past, each charge group was represented by the near- 
est mass on the stability curve for the associated frag- 
ment charge. The most recent versions of the trans- 
port computer codes use an isobaric flux representa- 
tion with the nearest charge on the stability curve; 
the distance to the nearest isobar was calculated as 

D = (Ai- A e ) 2 + 4 (Z { - Z ( ) 2 (44) 


where Z\ is the fragment and A ^ Z( is the listed 
isobar mass and nearest charge to the stability curve 
used in the calculation. In the following section, a 
brief description is presented of the nuclear models 
used in NUCFRAG. 

Total Absorption Cross Sections 

The nucleon-nucleus absorption cross sections are 
given by Letaw’s formula (ref. 22) that was con- 
structed to fit a consistent set of measurements made 
by Bobchenko et al. (See ref. 23.) For nucleus- 
nucleus collisions, an energy-dependent parameter- 
ization of a Bradt and Peters (ref. 24) form given by 
Townsend and Wilson (ref. 25) is used. Recently, a 
more accurate treatment for light ion-nucleon colli- 
sions has been derived in reference 26 in which a nor- 
malization factor of 0.95 is added to the parameter- 
ized expression of Townsend and Wilson at energies 
above 80 MeV/A. At energies below 80 MeV/A, a 
separate expression resulted from the complex quan- 
tum mechanical calculation given in reference 26. 

Fragmentation Cross Sections 

The fragmentation cross section of heavy ions 
is the least known physics quantity as inputs to 
the heavy ion transport calculations because of 


the dearth of experimental measurements. In the 
semiempirical nuclear fragmentation model devel- 
oped for NUCFRAG (ref. 27), the geometric abrasion- 
ablation model of Bowman, Swiatechi, and Tsang 
(ref. 28) was modified for the effect of frictional 
spectator interaction (FSI) by using a semiempirical 
correction to the abraded prefragment excitation 
energies. The important effect of the interacting 
electromagnetic fields (refs. 29 and 30) for heavy 
ions is also added. For nucleon-nucleus collisions, 
the fragmentation cross section is generated accord- 
ing to Silberberg and Tsao (refs. 31 and 32) rather 
than by the abrasion-ablation model. An improve- 
ment to fragmentation cross section for the light ions 
is given in reference 26. 

Renormalization 

The total absorption and fragmentation cross sec- 
tions are generated by different models, so a re- 
normalization process is needed to conserve mass 
and charge. Note that the earlier models of these 
cross sections failed to conserve mass and charge 
(ref. 33) and exhibited mass loss of up to 30 per- 
cent for 10 < Z < 22. The mass loss is displayed in 
figure 3 in which cr abs is compared with £ A,cr, p /A p , 
where A z is the fragment mass, o\ v is the fragmen- 
tation cross section for the projectile p, and A v is 
the projectile mass. Because the fragmentation cross 
section is less certain than the total absorption cross 
section, the former is renormalized to agree with the 
latter. 


Results 

The measured LET distributions behind 2 g/oii^ 
aluminum, 5-cm polyethylene, and 8-cm polyethylene 
targets arc shown in figures 4(a) 4(c), respectively. 
Also shown are the analytical results obtained by 
mapping the calculated LET distributions of attenu- 
ated iron projectile and projectile fragments into the 
detector response function where the premapping re- 
sults were obtained from a two-layer Green s function 
computer code (GRNTRN) assuming a 2.24 g/cm 2 
lead foil and target combination. Although the shape 
of the response function is not exactly known, a cor- 
rection for non-Gaussian contributions is taken as 


R(L, Lo) = 0.8 


+ 0.2 
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\/ 27nT o 
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exp 
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where a 0 = 0.4247r and cr\ = 2.4cro are fitted to 
the high-LET side of the measured primary ion peak 
for a 2 g/cm 2 aluminum target. (See fig. 4(a).) 
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Reasonable agreement between the theory and the 
measured LET distribution was obtained for all the 
targets; in general, the predicted iron peaks occur at 
slightly lower LET than the measured peaks, which 
indicates that the inferred energy of 557 MeV /A may 
be slightly high. Also, the predicted production of 
charge 24 and 25 ion groups are much higher than 
that of the other fragments. This discrepancy may 
have resulted from the choice of the current isotope 
set. (See table II.) For the overall fragment flux level, 
better agreement is noted for the polyethylene targets 
than for the aluminum target. The assumption that 
the fragmentation cross sections are too small for 
aluminum is inconclusive because the low production 
of fragments by 2 g/cm 2 aluminum can be more 
significantly affected by the fragment contribution 
originating from the media in front of the target 
than that currently considered in the analysis. For 
the high-LET side of the primary peak, the least 
agreement is noted for the 8-cm polyethylene because 
the non-Gaussian shape corrections were based on 
the aluminum target and the shape may vary for 
different LET regions. (Note that the peak for the 
thicker target appears at a higher LET.) 


Concluding Remarks 

A preliminary analysis of the measured LET spec- 
tra of a fragmented iron beam at 600 MeV/A nominal 
extraction energy at BEVALAC has been performed 
for 2 g/cm 2 aluminum, 5-cm polyethylene, and 8-cm 
polyethylene targets. Reasonable agreement is ob- 
tained between the experimental data and the cal- 
culations by using a two-layer (lead foil and target) 
Green’s function computer code. The comparison 
with experiment could possibly be improved by us- 
ing a more accurate isotope set or by including in the 
transport calculation the effect of other important 
beam-degrading media in front of the target. Fu- 
ture improvements in the code could include energy- 
dependent nuclear cross sections, the secondary frag- 
ment spectra, straggling, and multiple scattering but 
will have little effect on the results presented in this 
report. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
February 17, 1994 
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Table I. Summary of April 1992 Bevalac Irradiation 
600-Mev/A Iron at Extraction 



Density = 0.92 g/cm 2 . 


Table II. Isotopes List 
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Projectile ion charge 


Figure 3. Absorption cross sections in hydrogen component of polyethylene target and mass-average 
cross sections for various projectiles at 600 MeV/A. 
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Transported flux/incident flux 



LET, keV/pm 

(a) 2 g/cm 2 aluminum target. 

Figure 4. Comparison of two-layer (2.24 g/cm 2 lead and target) Green’s function calculation (GRNTRN) for 
557-MeV /A iron beam and measured LET distribution. 
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